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ABSTRACT 

The non-Gaussian cold spot found in the WMAP data has created controversy 
about its origin. Here we calculate the Bayesian posterior probability ratios for three 
different models that could explain the cold spot. A recent work claimed that the 
Spot could be caused by a cosmic texture, while other papers suggest that it could 
be due to the gravitational effect produced by an anomalously large void. Also the 
Sunyaev-Zeldovich effect caused by a cluster is taken into account as a possible origin. 
We perform a template fitting on a 20° radius patch centered at Galactic coordi- 
nates (b = —57°, I = 209°) and calculate the posterior probability ratios for the void 
and Sunyaev-Zeldovich models, comparing the results to those obtained with texture. 
Taking realistic priors for the parameters, the texture interpretation is favored, while 
the void and Sunyaev-Zeldovich hypotheses are discarded. The temperature decre- 
ment produced by voids or clusters is negligible considering realistic values for the 
parameters. 
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1 INTRODUCTION 

The study of the temperature anisotropics in the Cosmic 
Microwave Background (CMB) can provide crucial informa- 
tion on the origin of the Universe and help us to discriminate 
between different cosmological models. Recent observations 
' in many cosmological fields tend to support the standard 
inflationary model, which predicts the anisotropics to rep- 
resent an isotropic Gaussian random field. The Wilkinson 
Microwave Anisotropy Probe (WMAP) 1-year data (Ben- 
nett et al. 2003), measured the CMB fluctuations with high 
accuracy. In a first approach (Komatsu et al. 2003), these 
were found to be consistent with the Gaussian predictions. 
Further analyses revealed asymmetries or non-Gaussian fea- 
tures which have been confirmed in the 3-year WMAP data 
(Hinshaw et al. 2007) . Some of these anomalies are: low mul- 
tipole alignments (de Oliveira-Costa et al. 2006, Copi et al. 
2007, Land & Magueijo 2007, Chiang et al. 2007), north- 
south asymmetries (Eriksen et al. 2007), structure alignment 
(Vielva et al. 2007), non-Gaussian features found with steer- 
able (Wiaux et al. 2008) and directional wavelets (McEwen 
et al. 2006), and a very cold spot in the southern hemi- 
sphere (Cruz et al. 2007a). This cold spot (CS) is found 
at Galactic coordinates (b — — 57°, I = 209°), and has an 
angular radius in the sky of ~ 5°. The Spherical Mexican 
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Hat Wavelet (SMHW) was used to detect this structure. 
This spherical wavelet is an optimal tool to enhance non- 
Gaussian features on the sphere. The first SMHW analysis 
(Vielva et al. 2004) was performed on 15 different wavelet 
scales and two estimators, skewness and kurtosis, were used. 
The kurtosis deviated significantly from the range of val- 
ues expected from Gaussian simulations at scales around 
5° . The spot was found to be the main cause of the devia- 
tion from Gaussianity, using the area estimator (Cruz et al. 
2005). Cruz et al. (2006) showed the incompatibility of the 
frequency independent signal of the CS with Galactic fore- 
grounds, and Cayon, Jin & Treaster (2005) introduced two 
new estimators, namely Max and Higher Criticism, rede- 
tecting the CS. Cruz et al. (2007a) confirmed the detection 
in the 3-year WMAP data, using all the mentioned estima- 
tors. In the latter paper, the probability of finding an at 
least as high deviation in Gaussian simulations was calcu- 
lated without any a posteriori assumption. A conservative 
estimate of this probability was 1.85%, and hence this spot 
is a very special feature whose origin deserves further in- 
vestigation. Rath et al.(2007) detect the CS in the 3-year 
WMAP data using scalar indices instead of wavelets. 

Several models were proposed in order to explain this 
anomaly: voids (Inoue & Silk 2006, 2007, Rudnick, Brown & 
Williams 2007), second order gravitational effects (Tomita 
2005, Tomita & Inoue 2007), anisotropic Bianchi Vllh model 
(Jaffe et al. 2005), finite cosmology model, (Adler, Bjorken & 
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Overduin 2006), asymptotically flat Lemaitre-Tolman-Bondi 
model (Garcia-Bellido & Haugbolle 2008) and brane-world 
model (Cembranos et al. 2008) are some of them. Also a 
large-scale, non-Gaussian angular modulation (Naselsky et 
al. 2007) has been suggested although this hypothesis is not 
based on any physical model. The Bianchi VHh model was 
proven by Jaffe et al. (2006) to be ruled out at high sig- 
nificance, whereas there is still no further evidence for the 
validity of the other explanations. A full Bayesian evidence 
analysis of that Bianchi model can be found in Bridges et al. 
(2007a). Bridges et al. (2007b) show that the cold spot was 
likely to be driving any Bianchi Vllh template detection. 

Recently, Cruz et al. (2007b) showed that the CS could 
be caused by a cosmic texture. The amplitude and scale 
of the spot were consistent with that interpretation and the 
kurtosis of the data was compatible with the Gaussian CMB 
plus texture model at all scales. The predicted number of 
spots at a scale of 5° or larger, is ~ 1 which is consistent 
with the observed spot. 

In this paper we study two other possibilities, namely 
the Rees-Sciama (Rees & Sciama 1968) effect caused by 
voids as proposed by Inoue & Silk (2006, 2007), Rudnick, 
Brown & Williams (2007) and the Sunyaev-Zeldovich effect 
produced by clusters (hereafter SZ effect, Sunyaev & Zel- 
dovich 1970) already considered in Cruz et al. (2005). Both 
of them could in principle produce temperature decrements 
as the CS. In addition a dip (Rudnick, Brown & Williams 
2007) in the NRAO VLA Sky Survey data (NVSS, Condon 
et al. 1998) and the Eridanus group of galaxies (Brough et al. 
2006) match approximately the size and position in the sky 
of the CS. However, the values required by the parameters 
of either model are not very realistic. A Bayesian framework 
is very convenient to decide which model best describes the 
data. The prior theoretical and observational knowledge on 
the parameters is included in the Bayesian analysis and can 
be decisive in some cases. Our analysis is performed on the 
recently released 5-year WMAP data (Hinshaw et al. 2008) . 
The large-scale structure of the CMB anisotropies remains 
almost unchanged respect to the 3-year release (Komatsu et 
al. 2008). 

We will first review the data and the template fitting 
method in section 2, then we recall in section 3 the results 
obtained for the texture model in order to compare with the 
SZ effect in section 4 and the void model in section 5. Finally 
we present our conclusions in section 6. 



2 BAYESIAN TEMPLATE FIT 

We use the combined, foreground-cleaned Q-V-W map 
(hereafter WCM) of the 5-year WMAP data release (Hin- 
shaw et al. 2008). The WCM is in the HEALPix pixelization 
scheme (Gorski et al. 2005) with resolution parameter Nside 
= 64. Since the scale of the spot is around 5° this resolution 
is sufficient and reduces the number of pixels used in the 
template fitting. 

The template fitting is performed in a circular area of 
20° radius centered at Galactic coordinates (b = — 57°, I = 
209°). Although the angular radius of the CS is about 5° 
we have to consider at least a 20° radius patch to take into 
account the whole neighborhood of the spot for the fit. The 



SMHW convolves all the pixels in this region and they con- 
tribute in an important way to the detected structure. 

In our template fitting we apply Bayesian statistics to 
the data, D, in order to find the best Hypothesis, H, de- 
scribing the data through a set of parameters, 0. 

Bayes' theorem states: 



P(&\D,H) = 



P(D\&,H)P(&\H) 
P(D\H) 



(1) 



where P(®\D,H) is the posterior distribution, P(D\&, H) 
the likelihood, P(&\H) the prior and P(D\H) the Bayesian 
evidence. 

The null hypothesis, Ho, describes the data as an 
isotropic Gaussian random field (CMB) plus the noise of 
the WMAP data. The alternative hypothesis, Hi describes 
the data as CMB plus noise and an additional template, T, 
given by a physical model as for instance a texture spot, the 
SZ effect produced by a cluster or the Rees-Sciama effect 
generated by a local void. Hi (with i = 0,1) depends on the 
parameter set &i. 

For Hi, the probability density associated with the ob- 
served data is 



P(D\Hi) 



P(D\Bi,Hi)P(®i\Hi)dei. 



(2) 



The evidence is hence the average of the likelihood with re- 
spect to the prior and has been extensively used for model 
selection in cosmology (e.g. Jaffe 1996, Hobson, Bridle & La- 
hav 2002, Hobson & McLachlan 2003, Mukherjee, Parkin- 
son & Liddle 2005, Parkinson, Mukherjee & Liddle 2006). 
A model having a large evidence shows a high likelihood for 
its allowed parameter space given the data. 

We redefine our notation, with the likelihood L, and the 
prior n so the previous integral for the evidence, E, reads: 



E t = / L(&i\Hi)U(&i)d&i. 



(3) 



The likelihood function is L oc exp (— x 2 /2), where \ 2 
is defined as: 



X 



(D — T) N (D — T), 



(4) 



and N is the generalised noise matrix including CMB and 
noise. The template T depends on the parameter set ©. The 
calculation of the noise contribution to N is straightforward 
since the noise in the WCM is uncorrelated and well known. 
In order to obtain the CMB part, we calculate the corre- 
lation function for the WCM taking into account the pixel 
and beam effects. Comparing the resulting evidence ratios 
obtained from the analytical expression for the WCM cor- 
relation function with those obtained using simulations, the 
errors are negligible. 

The posterior probability ratio: 



Pr(gi|D) = Ei Pr(ffx) 
Pr(# |D) E Pt(Ho)' 



(5) 



can be used to decide beween both hypotheses. The alter- 
native hypothesis is the more probable one when p > 1. If 
p < 1 the null hypothesis is the preferred one to describe 
the data. The a priori probability ratio for the two mod- 
els, Pr(Hi)/ Pr(Ho) is set to the fraction of the sky, fs, 
covered by either texture, clusters or voids. In this man- 
ner we compensate the a posteriori selection of the pixel 
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where the template is centred, namely the center of the CS, 
(b = -57°, I = 209°). 



3 TEXTURE 

Unified theories of high energy physics predict the produc- 
tion of topological defects (Kibble 1976, Vilenkin & Shellard 
2000) during a symmetry-breaking phase transition in the 
early universe. Depending on the order of the broken sym- 
metry, different types of defects arise, such as domain walls, 
strings, monopoles or texture. Each of these extremely ener- 
getic events leave a characteristic pattern in the CMB. For 
instance cosmic texture (Turok 1989) produce hot and cold 
spots in the CMB sky. 

Cruz et al. (2007b) showed recently that the CS could 
be due to a collapsing texture. A Bayesian template fit was 
performed on the 3-year WMAP data, using an approxi- 
mated temperature profile based on the analytic result given 
in Turok & Spergel (1990). The a priori probability ratio for 
the two models, Pr(i3i)/ Pr(ffo), was estimated in this case 
by the fraction of sky covered by texture, fs = 0.017. The 
evidence ratio was Ei/Eo ~ 150, hence the resulting poste- 
rior probability ratio was p = 2.5 favoring the texture hy- 
pothesis versus the null hypothesis. The scale and amplitude 
of the spot (5.1° and dT/T = 7.7 x 10~ 5 ) were consistent 
with having being caused by a texture at redshift z ~ 6 and 
energy scale <f>o ~ 8.7 x W 15 GeV. The number of texture 
spots of scale 5.1° or above is around 1.1 consistent with 
the single observed spot. 

The amplitude of the texture spot was shown to be over- 
estimated by a selection bias. Texture spots with amplitude 

4 x 10~ J were simulated and added to CMB Gaussian sim- 
ulations. Performing a template fit on these single extreme 
events, the inferred amplitude was in average 7.9 x 10~ 5 , i.e. 
almost twice the true amplitude. This means that the am- 
plitude required by the CS is not in conflict with the upper 
limits on the amplitude calculated from the power spectrum 
(Urrestilla et al. 2007). 

Moreover, Bevis et al. (2008) and Urrestilla et al. (2007) 
find a 2<j preference for a cosmological model including tex- 
ture or strings, when considering only CMB data. How- 
ever this preference is reduced when including Hubble Key 
Project (Freedman et al. 2001) and Big Bang nuclesosynthe- 
sis (Kirkman et al. 2003) constraints in the Bayesian power 
spectrum anlysis. 

In order to estimate the significance of the p = 2.5 re- 
sult, the same fitting procedure was performed on 10,000 
Gaussian CMB simulations. The most prominent spot was 
selected for each simulation and a texture profile was fitted 
to it. Around 5.8% of the simulations showed p > 2.5 and 
the typical values where p « 0.14. This result was found to 
be almost independent on the priors or the exact shape of 
the approximated temperature profile. 

Using the 5-year data release we obtain Ei/Eo « 160 
and p — 2.7. The lower noise level allows a slightly better 
detection of the texture template in the new data release. 
We have not repeated the frequentist estimation of the sig- 
nificance because the result was very similar to the 3-year 
one. The probability of having an at least as extreme value 
(i.e. p- value) of p assuming that the null hypothesis is true, 
should be again around 5%. 



The presence of texture would explain the excess of kur- 
tosis observed at scales around 5°. The 1° — 2° texture pre- 
dicted by the model would be confused with Gaussian CMB 
spots since the latter have their maximum power at these 
scales. 

A number of follow-up tests can be performed to check 
the texture hypothesis. For instance an all-sky template fit- 
ting looking for several texture can confirm or refute the 
alternative hypothesis. Also a polarization analysis can help 
to decide, since Gaussian CMB spots have a correlated po- 
larization pattern, whereas late texture do not show asso- 
ciated polarization. In addition the B-mode texture power 
spectrum is expected to be larger than the typical inflation- 
ary one, having also a different shape (Urrestilla et al. 2007). 



4 CLUSTERS 

The hot gas trapped in clusters modifies the CMB spectrum 
through the Sunyaev-Zeldovich effect. The high-energetic 
electrons present in the gas collide with the CMB photons, 
changing their original frequency. Two different components 
can be identified. The former is produced by random ther- 
mal motion of the electrons (thermal effect), and the second 
one, by the peculiar motion of the cluster with respect to 
the CMB (kinetic effect). The thermal SZ effect is typically 
one order of magnitude higher than the kinetic one, and 
therefore we focus on the first one. The effect on the CMB 
photons is frequency dependent and at frequencies below 
217 GHz it produces a temperature decrement. Hence this 
effect is a possible candidate to explain an anomalously cold 
spot at the WMAP frequencies which range from 23 GHz to 
94 GHz. 

The temperature decrement is given by: 

AT 

-jr=g{x)y c , (6) 

where x = hv/kBTo is the frequency in dimensionless units, 
with h and fes Planck and Boltzmann constants respectively, 
and g(x) = (xcoth(a;/2) — 4). The cluster Comptonization 
parameter is defined as 

vrks r Te (7) 

m e c 2 J 

where n e and T e are the intracluster electron density and 
temperature, ot is the Thomson cross section, and m e is 
the electron mass. The integral is performed along the line 
of sight through the cluster. In hot clusters (T e > 5 keV), 
the relativistic electrons lead to a slightly modified spectral 
behaviour of the thermal SZ effect, but this is negligible in 
our analysis. Also neglected is the kinetic SZ effect of the 
cluster. 

Cruz et al. (2006) showed that the flat frequency de- 
pendence of the CS is incompatible with being caused by 
a SZ signal alone. However a combination of CMB plus SZ 
effect may explain the spot and could have a sufficiently flat 
frequency spectrum. 

Cruz et al. (2005) already reported that a group of lo- 
cal galaxies were present in the direction of the CS. They 
argued that the amount of gas needed to explain the whole 
signal was too large for a small group of galaxies. However 
the group of galaxies in the direction of the spot has been 
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Figure 1. Azimuthal projection of an pa 18° X 18° patch of the 
WMAP combined map (WCM, top) and the Nearby Galaxies 
Catalog (bottom), centered at (b = —57°, I = 209°) and con- 
volved with the SMHW at scale 5° . The units are wavelet coeffi- 
cients of temperature in mK and number of counts for the WCM 
and Nearby Galaxies Catalog respectively. 



shown recently to be larger than previously thought. Brough 
et al. (2006), define it as the Eridanus supergroup with mass 
M ta 1O 14 M0. The angular size and position in the sky of 
Eridanus are similar to those of the CS, as can be seen in 
Figure (1). The top panel shows the WCM data and the bot- 
tom one, the Nearby Galaxies Catalog (Tully, 1987). Both 
are as 18° x 18° patches centered at (b = -57°,/ = 209°) 
and were convolved with the SMHW at scale 5°. However 
there are other, more prominent nearby groups which do not 
have an associated CMB spot, hence this could be a statis- 
tical coincidence, but it is worth testing through a template 
fit whether there is a hidden SZ-like template in the WCM 
data. 

The temperature profile as a function of the angle in 
the sky, 6, can be calculated assuming a /3-model (with f3 = 



3/2) for the gas-density profile (see e.g. Atrio-Barandela & 
Miicket, 1999): 



Ar) = n c (l + (r/r c ) 2 )-^ /2 , 



(8) 



where n c is the central electron density and r c is the core 
radius of the cluster. Expressing the virial radius of the clus- 
ter r v in terms of the core radius as r v = pr c , Eq. [7] can be 
rewritten as: 



where y = {Hbctt /m e c 2 ]r c Tn c , and 



V 1 + (W 



: tan 



(0/0c 



i + (e/e c 



(9) 



(10) 



The parameter 6 C is the angle subtended by the core radius 
of the cluster, whereas p is not a parameter but a function 
of mass M and redshift z. Typical values for clusters yield 
p = 10 - 15. 

Summarising, our temperature profile with parameters 
f? c and y reads: 



AT 



tan 



P 



(W 2 



i + (o/e c y 



(ii) 



Eq. [11] gives the effect produced by a single cluster through 
a particular line of sight. As the WCM is a combination 
of maps at different frequencies, we choose the central fre- 
quency, v — 61 GHz. 

We consider two alternative hypotheses in this context. 
In the first one the spot would be caused by a rich cluster 
and in the second one by a nearby group of the Nearby 
Galaxies Catalog. 

We simulate the entire population of clusters with flux 
above 10 mjy (at 353 GHz), using a Monte Carlo algo- 
rithm where masses and redshifts were assigned to each 
cluster according to the Press-Schechter prescription (Press 
& Schechter 1974). We assume Q. m = 0.3, Oa = 0.7 and 
as = 0.8 (where Q m and Q,a are the matter and vacuum 
energy density parameters, and as the dispersion of the den- 
sity field at scales of 8 Mpc) . The cut in flux removes nearby 
small groups of galaxies and distant small clusters. The red- 
shifts of the simulated clusters range from 0.005 to 3.4. For 
each cluster with known mass and redshift the core radius 
is determined from the scaling relation: 

r c = 150Mj/ 3 (l + 2)" 1 fe~ 1 kpc (12) 

where Mis is the mass of the cluster in 10 15 /i _1 Mq. It is 
important to note that the cut in flux does not have any neg- 
ative effect in the results of this paper since the fluxes we are 
trying to explain are much larger than 10 mjy. Since the red- 
shift and core radius are known for each cluster, we calculate 
the angle subtended by the core, C . We find the maximum 
angle of the 300,000 clusters to be 6 C » 0.2°. Therefore the 
cluster hypothesis has to be discarded for generating a 5° 
spot. 

Returning to the particular case of the Eridanus super- 
group found in the Nearby Galaxies Catalog, we can consider 
larger values of 8 C , due to its proximity and lower values of 
y since it is not a cluster but a big group or an ensemble 
of groups. In addition we assume lower values of p because 
the gas could be rather dispersed since the gravitational at- 
traction is lower than in big clusters. Now the alternative 
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hypothesis is that a big group from the Nearby Galaxies 
Catalog is causing the CS. 

Although the center of Eridanus does not match ex- 
actly the center of the spot, we still center the template at 
(6 = —57°, I = 209°), since the gas does not need to fit ex- 
actly the positions of the galaxies. This assumption is the 
most favorable for the alternative hypothesis. If the CS was 
due to the SZ effect we would expect high X-ray emission in 
the considered region. We can use the ROSAT all-sky sur- 
vey (Snowden et al. 1997) to estimate the prior on y . The 
ROSAT survey excludes some small regions of the sky, and 
one of them crosses the considered patch. However the es- 
timated average counts in the ROSAT R6 band are around 
80 cts s _1 arcmin -2 . The /3-models which are compatible 
with this background and at the same time show an ex- 
tended SZ signal of several degrees, range between models 
with a core radius of 0.05 Mpc and bolometric luminosities 
of 10 42 erg/s and models with 0.15 Mpc and bolometric lu- 
minosities of 7 x 10 42 erg/s. These models would predict at 
most y = 10 -7 . The nearest groups in the Nearby Galaxies 
Catalog are at 5 Mpc distance and as we have shown, the 
upper limit for the core radius is 0.15 Mpc. In this case we 
get 8 C = 0.03 rad, which we set as the upper limit for 8 C . 
We set the lower limit on 6 C to 1°, since the spot has 5° 
radius and the pixel size is around 1° for the chosen res- 
olution. Hence we consider flat priors for parameter y in 
the interval [0, 10" 7 ] and for 6 C in [0.0175, 0.0300] radians. 
Setting p = 5 we obtain a relative evidence of Ei/Eq ~ 1. 
The fraction of the sky covered by nearby galaxies in the 
catalog, is fs ~ 0.05 therefore p ta 0.05 < 1 and the SZ 
hypothesis is discarded again. The best fit parameters are 
the upper limits of the priors, y « 10 -7 and 8 C = 0.03 
rad. The amplitude of the temperature decrement does not 
exceed lfiK and is hence negligible compared to the deep 
decrement. Therefore the E\/Eq « 1 result is not suggesting 
that the SZ hypothesis is viable for explaining the spot. As 
the template has a negligible effect on the data, T « 0, the 
likelihood is almost constant in the whole parameter range 
and as the priors integrate to unity we have E\/Eq ~ 1. 
This is a well-known feature of Bayesian model selection. 
In any problem, if H\ has additional parameters compared 
to Ho, and the likelihood is insensitive to these parameters, 
then Ei — Eq. The results do not significantly depend on 
the value of p, the chosen priors or the resolution. 

The gravitational effect due to a supercluster cannot 
strongly affect large angular fluctuations because these mas- 
sive structures are very concentrated in a particular region. 
On the contrary, voids are better candidates to generate 
large angular spots. 



5 VOIDS 

The gravitational effect produced by cosmic voids can pro- 
duce temperature decrements in the CMB and has been pro- 
posed as a possible source of the CS (Inoue & Silk 2006, 
2007). In analogy to the cluster case with the Eridanus su- 
pergroup, there is a dip in the NVSS radio source catalogue 
that approximately matches the position of the spot. Rud- 
nick, Brown & Williams (2007), claim that the dip could be 
associated with a void at redshift z < 1 which could be the 
origin of the CS, but affirm that the significance of this co 
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Figure 2. Azimuthal projection of an pa 18° X 18° patch of the 
WCM (left column) and the NVSS (right column), centered at 
(b = -57°, I = 209°) and convolved with the SMHW at scale 
1.7° (top row) and scale 5° (bottom row). The units are wavelet 
coefficients of temperature in mK and radio source number counts 
for the WCM and NVSS respectively. 



incidence is unclear. As the fraction of the Universe filled by 
voids is around 30% in volume at redshift z fa 1 (Colberg et 
al. 2005), the fraction of the sky covered in projection, fs, is 
close to one and hence such coincidences are not rare. The 
dip and the CS appear at different scales as can be seen in 
Figure 2 where we present a « 18° x 18° patch of the WMAP 
and NVSS data convolved with the SMHW at scales 5° and 
1.7°. At scale 5° the CS reaches its maximum power in the 
WCM, whereas there is no significant dip visible in the NVSS 
data. At scale 1.7° where the CS is not an outstanding fea- 
ture, a dip becomes visible but there are other similar dips 
with no associated CMB spots, therefore the association be- 
tween dip and spot is unclear and could be a coincidence. 
Moreover, there are several other dips with similar depth 
outside the selected region hence the dip does not seem to 
be an outstanding feature of the NVSS data. McEwen et 
al. (2007) found several regions where both data sets were 
highly correlated and one of them was located close to the 
CS. This could be due to the very cold temperature of the CS 
which gives a high correlation with any underdense region in 
the NVSS data at scale 5° and not to the relation between 
dip and spot since they have different scales (see Figure 2). 
Rudnick, Brown & Williams (2007) also performed an order 
of magnitude calculation to estimate the radius of the void, 
r v , needed to generate such a CMB spot through the Rees- 
Sciama effect obtaining r v ~ 200/i _1 Mpc in concordance 
with previous calculations (Inoue & Silk 2007). The prob- 
ability of observing such a void in standard theory is very 
small. For instance, a void with r v = 140/i _1 Mpc will occur 
with a probability of 5 x 10~ 10 in an optimistic case (Rud- 
nick, Brown & Williams 2007). Typical observed values are 
r v ~ 10/i -1 Mpc for a matter density contrast of S m ~ —0.8 
(Hoyle & Vogeley 2004, Patiri et al. 2006). Similar results 
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are found in N-body simulations (Colberg et al. 2005). How- 
ever it is worth performing a Bayesian template fit for the 
void hypothesis. 

Martinez-Gonzalez et al. (1990) and Martinez- 
Gonzalez & Sanz (1990), performed a calculation of the tem- 
perature anisotropy generated by a compensated void in an 
Einstein-de Sitter Universe, 



Prior on rv 



AT _ 16 rj 
T ~ 3 (1-dcosa) 3 
9 



x 



- 



d X .2 
sin a 



(13) 



where a is the angle formed between the direction of obser- 
vation and the direction defined by the centre of the struc- 
ture, and d is the distance to the centre of the structure. 
They adopt the thin shell approximation, where the void is 
an empty spherical region surrounded by a thin shell con- 
taining the matter which would have been inside the void. 
Outside the void the Universe is represented by an Einstein- 
de Sitter model and the shell propagates approximately as 
r v oc f and 7 ~ 0.8. The units are 8ttG = c = dh = 1, 
where dh is the horizon distance. Eq. [13] applies to photons 
crossing the void. The photons that do not cross the void 
are unaffected since there is no net mass excess. 

It is straightforward to show that Eq. [13] is still valid in 
a flat, A - Cold Dark Matter (A-CDM) Universe, assuming 
that r v oat 1 , and considering that 

dz (14) 

H F( z y ll4j 
F(z) = y/n A + (1 + z ym m + (i + z ) 4 n r , (15) 

where Q, r — 10 -4 is the radiation density parameter. If we 
further assume that r v remains constant in comoving coor- 
dinates we have: 

dz 



d = 



7 



d(log a) 



F(z) 



(1 



(16) 



d(logt) 

with scale factor a. 

Introducing the above expressions for d and 7 in Eq. [13] 
we have our temperature template with parameters z, r v in 
a flat A-CDM Universe. 

We choose the prior on r v in concordance with the void 
radius distribution found in the Two-Degree Field Galaxy 
Redshift Survey data (Fig. 3 in Hoyle & Vogeley 2004). We 
extrapolate these values for large r v through an exponential 
decay as shown in Figure (3). This extrapolation is gener- 
ous with big voids since the derived decay is faster than the 
exponential used. For the redshift, we choose a flat prior 
between and 1. This assumption is again rather generous 
with large radii, since according to N-body simulations (Col- 
berg et al. 2005) the void sizes should be smaller at high, 
unobserved redshifts. We set the a priori probability ratio 
for the two models to unity since voids cover almost all the 
sky. The obtained evidence ratio is E1/E0 ~ 1 and therefore 
the posterior probability ratio is p ~ 1. This result is not 
suggesting that the void model is viable as the best fit tem- 
plate corresponds to r v = 27h~~ 1 Mpc and z = 0.03 which 
produces a temperature anisotropy of 0.2pK. As in the clus- 
ter case the template is completely negligible, the likelihood 
almost constant, and the priors integrate to unity, therefore 
E1/E0 ~ 1. The contribution of voids with r v > 30/i -1 Mpc 




Radius of the void h , Mpc 

Figure 3. Prior on the radius of the void parameter (red line). 
The underlying histogram is taken from the distribution found 
in the Two-Degree Field Galaxy Redshift Survey data (Fig. 3 in 
Hoyle & Vogeley 2004). 



to the evidence integral is negligible and the results remain 
unchanged varying the upper limit of parameter r v for val- 
ues above 30/i _1 Mpc. 

The big voids (r v > 100/i~ 1 Mpc and z < 1) described 
in Rudnick, Brown & Williams (2007) which do produce an 
important temperature decrement are strongly penalized by 
the prior, despite showing higher likelihood values. We re- 
define Hi in order to test the big void hypotesis. Now only 
voids with r v > 100/t _1 Mpc and z < 1.0 are considered. 
We renormalise the priors in the new parameter ranges and 
recalculate the prior probability ratio since the big voids 
(in case they exist) are expected to cover only a small 
part of the sky. As shown above a void with r v — 140h~ 
Mpc is expected to occur with a probability of 5 x 10~ 10 
(Rudnick, Brown & Williams 2007). Hence as a valid up- 
per limit we can assume that there is at most one big void 
with r v — 140/i _1 Mpc and 2 = 1. The fraction of the sky 
covered by this void is fs ~ 0.008. Therefore we assume 
Pr(fli)/Pr(#o) = 0.008. Recalculating the evidence for the 
new parameter range, we obtain E1/E0 « 2 which yields 
p ~ 0.016. This result discards the possibility of explaining 
the CS through such a void. Hence the void hypothesis has 
to be rejected. 

Inoue & Silk (2006, 2007) proposed slightly different 
voids with very high values for the radius of the void, r v ~ 
300/i _1 Mpc and quasi linear density contrast values 5 m w 
-0.3, while Rudnick, Brown & Williams (2007) used 5 m ~ 
— 1. Inoue & Silk (2007) recognize that the required values 
for r v are inconsistent with observations at low z and N- 
body simulations. Such a void is a 13cr result in the standard 
theory with as = 0.74. They consider that a percolation 
process of the typically observed small voids could form the 
large quasi linear voids. However these voids have not been 
observed yet and as already mentioned, according to N-body 
simulations (Colberg et al. 2005) the void sizes should be 
smaller at high, unobserved redshifts. 
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6 CONCLUSIONS 

In this paper we consider different possibilities to explain a 
non-Gaussian cold spot found in the WMAP data. This spot 
is located at Galactic coordinates (b = —57°,/ = 209°), cov- 
ers around 10° in the sky, has been shown to be incompatible 
with Gaussian simulations and is unlikely to be produced by 
foregrounds (Vielva et al. 2004, Cruz et al. 2005, Cruz et al. 
2006, Cruz et al. 2007a). We compare three alternative hy- 
potheses to explain the CS, namely texture, clusters and 
voids versus the null hypothesis. The latter describes the 
data as an isotropic Gaussian random field (CMB) plus the 
noise of the WMAP data. A Bayesian analysis takes into ac- 
count the prior observational and theoretical knowledge on 
the parameters and can be used to find the hypothesis best 
describing the data. The posterior probability ratio, p, is the 
quantity we use to compare the different hypotheses. The 
alternative hypothesis is preferred if p > 1 and disfavored 
otherwise. Cruz et al. (2007b) performed such an analysis for 
the texture model, finding that this model is favored versus 
the hypothesis of Gaussian fluctuations with p « 2.5. The 
predicted size, amplitude, temperature profile, and number 
of spots generated by texture is consistent with the observed 
spot. Repeating this analysis for the 5-year WMAP data we 
find p ~ 2.7 due to the lower noise level. The frequentist 
p-value associated to this result is around 5%. 

The finding of a local supergroup (Brough et al. 2006) 
and a dip in the NVSS survey (Rudnick, Brown & Williams 
2007), which match the position of the CS could be a hint 
for the cluster or void hypotheses, but our detailed Bayesian 
analysis reveals that these are statistical coincidences. The 
fraction of the sky covered by clusters and voids is included 
in our analysis. 

The large angular size and temperature decrement of 
the CS are hard to explain with conventional clusters or 
voids. We would need either a local, rich cluster or a void 
with radius r v ~ 100— 200/i _1 Mpc located at redshift z < 1. 
In the cluster case only a local group with Compton pa- 
rameter y < 10~ 7 is observed in the required direction. 
The angular scale is marginally consistent with the spot 
but the temperature decrement generated is only a few fiK 
and hence cannot account for the CS. Our Bayesian analysis 
gives p 0.05, discarding the cluster hypothesis. 

For conventional voids with r v ~ 10/i -1 Mpc we ob- 
tain p « 1 because the likelihood is insensitive to the extra 
parameters. A detailed analysis reveals that these voids pro- 
duce negligible temperature decrements below lpK. Renor- 
malising the priors for voids with radius r v ~ 100 — 200/i _1 
Mpc, we obtain a posterior probability ratio of p = 0.016. 
Therefore the void hypothesis is not viable either. Explain- 
ing the cold spot by a void with r v > 100ft -1 Mpc, would 
introduce a 13a anomaly in the standard model to solve the 
ft* 99% anomaly of the spot. The solution would be worse 
than the problem. 

Hence we can conclude that the texture hypothesis is 
the most plausible explanation for the CS. Many scientists 
are uncomfortable with cosmic defects, because these failed 
to explain structure formation. However, cosmic defects are 
a generic consequence of symmetry breaking, which is an 
essential ingredient in unified theories of particle physics. 
Hence it is very plausible that they exist. They might be 
formed at any scale below the Grand Unified scale, and 



would be consistent with current cosmological observations. 
They are even compatible with inflation. The only constraint 
is that they should make a subdominant contribution to 
structure formation (see for instance Urrestilla et al. 2007). 

The virtue of the texture hypothesis is that it is fully 
testable. A texture would not lead to a correlated polariza- 
tion pattern around the cold spot, but it would lead to a 
very specific gravitational lensing pattern. In addition there 
should be more CMB texture spots at scales around 1 — 2°. 
Although they are more difficult to detect because the Gaus- 
sian CMB reaches its maximum power at these scales, we 
may be able to detect them through an all-sky Bayesian 
analysis, which is already in progress. 
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